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Abstract 

We extend earlier work on leading logarithms in the massive nonlinear 0{n) sigma 
model to the case of SU(N) x SU(N)/SU(N) which coincides with mesonic chiral 
perturbation theory for N flavours of light quarks. We discuss the leading logarithms 
for the mass and decay constant to six loops and for the vacuum expectation value 
(qq) to seven loops. For dynamical quantities the expressions grow extremely large 
much faster such that we only quote the leading logarithms to five loops for the 
vector and scalar form factor and for meson-meson scattering. The last quantity we 
consider is the vector-vector to meson-meson amplitude where we quote results up 
to four loops for a subset of quantities, in particular for the pion polarizabilities. As 
a side result we provide an elementary proof that the factors of N appearing at each 
loop order are odd or even depending on the order and the remaining traces over 
external flavours. 

Keywords: Renormalization group evolution of parameters; Spontaneous and radia- 
tive symmetry breaking; Chiral Lagrangians; Meson-meson interactions 



Contents 



1 Introduction 1 

2 iV-navour mesonic Chiral Perturbation Theory 2 

3 Leading logarithms 4 

4 Mass 6 

5 Decay constant 6 

6 Vacuum expectation value 8 

7 Vector form factor 10 

8 Scalar form factor 14 

9 Meson-meson scattering 16 

10 77 — > Ten and pion polarizabilities 22 

11 Conclusions 28 
A Powers of N in the results 29 



1 Introduction 

The calculation of higher loop corrections is an important problem in all areas of particle 
physics. The leading logarithms in a renormalizable field theory can be calculated to all 
orders by simply using the renormalization group. In nonrenormalizable effective field 
theories like Chiral Perturbation Theory (ChPT) [1-3], the recursive argument underlying 
the renormalization group does not work since one has a new Lagrangian at each order. 
Weinberg [1] showed that using the requirement that all nonlocal divergences cancel, one 
could obtain the leading logarithms (LL) at two-loop order with only one-loop calculations. 
This method has then been applied to various processes at the two-loop level [4]. That 
it works to all orders was later proven using beta-functions [5] and also with a more 
diagrammatic method [6]. 

Using this method, [7-10] found recursion relations valid in the massless limit and 
applied them to a number of processes. Away from the massless limit the tadpoles do 
not vanish and this causes the number of needed one-loop diagrams at every order to 
increase considerably. A systematic method to automatize the calculational process was 
found in [6] and then applied to a number of processes in the normal [6, 11] and abnormal 
or anomalous intrinsic parity [12] sector of the massive 0(n) nonlinear sigma model. In 
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the present paper we extend the calculations in the even sector to the symmetry breaking 
pattern of SU(N) x SU(N)/SU(N). All results are for the case of equal masses. 

We discuss the leading logarithm contribution to the mass, decay constant, and vacuum 
expectation value to sixth or seventh order. Numerical results are discussed for the two 
physical cases N = 2,3. For the vector and scalar form factors we give expressions for 
the full results and for the radius and curvature. We present no numerical results, but 
some discussion of numerics for the scalar form factor for N = 2 can be found in [11]. For 
meson-meson scattering we present analytic results for the amplitude and the scattering 
lengths up to fifth order. We show numerical results only for the singlet scattering length 
for N = 3, which we compare with the full two-loop calculation as well. For 77 — > tttt 
we give analytic results for the full amplitude for general N and for the polarizabilities for 
N = 2. For the latter we also present numerical results. 

We provide some references to the N = 2 and N = 3 cases at the two-loop level where 
the general- N case is not known to that level. An extensive discussion of the corresponding 
literature can be found in the review [13]. 

In Sect. 2 we present the model and the different parametrizations we use. Sect. 3 
describes the changes needed compared to the 0(n) work and provides the necessary defi- 
nitions such that the formulas in this paper are self-contained. We do however not discuss 
in detail the methods used. The remaining sections present results for the mass (Sect. 4), 
decay constant (Sect. 5), vacuum expectation value (Sect. 6), vector form factor (Sect. 7), 
scalar form factor (Sect. 8), meson-meson scattering (Sect. 9) and vector- vector to meson- 
meson scattering (Sect. 10). In addition we prove in Appendix A that only certain powers 
of N can show up at each order. 

2 A^-flavour mesonic Chiral Perturbation Theory 

The Lagrangian of the massive nonlinear SU(N) x SU(N)/SU(N) sigma model or N- 
flavour mesonic ChPT at lowest order is given by 



where U is a special unitary N x N matrix, which contains N 2 — 1 degrees of freedom. 
(A) = tr(A). The interaction with external axial- vector and vector fields enters through 
the covariant derivative 



while the explicitly chiral symmetry breaking terms as well as the scalar and pseudoscalar 
external sources are contained in 



C = —(D,UD^) + —{ X tf + U X ] ) , 




(1) 



D ll U = dpU-±[v lli U\-±{a„U} i 



(2) 



X = 2B(s + ip) + M 2 1 . 



(3) 
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The chiral SU(N) x SU(N) symmetry is broken spontaneously to SU(N) by the vacuum 
expectation value (0|£/|0) = 1, where 1 is the N x N unit matrix. This leads to the appear- 
ance of N 2 — 1 Goldstone bosons, which correspond to the degrees of freedom contained 
in the matrix field U. The term proportional to M 2 breaks the symmetry explicitly and 
causes the Goldstone bosons to pick up a mass which, at tree level, is equal to M. In terms 
of equal quark masses rh we have M 2 = 2 Bin. 

The Lagrangian (1) coincides with ChPT and therefore constitutes an effective La- 
grangian for two- and three-flavour QCD for N = 2 and N = 3, respectively. Note, 
however, that in the case considered here, all mesons have the same mass. How this corre- 
sponds to a theory formulated in terms of quarks can be found in more detail in, e.g., [14]. 
Below we occasionally use a vector notation for quarks q with q T = . . . , q N ) where the 
subscript denotes the flavour. 

In previous publications [6, 11, 12], the chiral logarithms of the massive nonlinear 
0(n + l)/0(n) model have been considered. The two models coincide for N = 2 and 
n — 3, such that the corresponding results can be used as a check. 

There are many ways the special unitary matrix U can be parametrized in terms of the 
meson matrix (f> = <ft a T a , where T a are the generators of SU(N) normalized as (T a T b ) = 5 ab . 
Physical results are independent of this choice. As in the earlier work on the massive 0(n) 
model, one can therefore use different parametrizations to obtain a thorough check of the 
calculation. The four parametrizations we have used are 

.iV2\ l + >0 + ^) 

U\ = exp ——(f) , U 2 = 



V2F' 



1 + 




Us = e ' ft T3^r > ^ = e ^ \\l 1 ~^ 2 + i2 £* I • (4) 

The matrices must be special, i.e., dett/j = 1 which for U\ is an automatic consequence 
of ((p) = 0. For the other cases one has to solve for in terms of 0. Using (log Ui) = 
log det Ui = 0, one finds that the start at order 3 . For /3 2 we can then solve the resulting 
equation perturbatively while (3^ and (3^ can be written explicitly in terms of <fi as 

00 I i\n / / jl \ 2n+l^ 



N \ b V V2FJ b V V2F J/ N 



v (-i) w ( v 

^2n + l\ \V2Fj 



2n+l\ 



^ = ~Jr ( arCSin \^¥) ) = -Jf g 4"(n!H2n + 1) ( [f) /' (5) 

Note that the n = term vanishes in both sums since ((f)) = 0. We could have added a fifth 
parametrization by adding a singlet component to (f> in as was done for parametriza- 
tion 2. 

It is also possible to treat U(N) x U(N)/U(N) by simply allowing to have a singlet 
component and removing the /3j. We do not discuss this case. 
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3 Leading logarithms 



The method used here is entirely analogous to the work in [6, 11, 12] but with a traceless 
N x N Hermitian matrix instead of an n-dimensional vector. The calculations are done 
schematically as follows: first we generate all needed one- loop diagrams with a CH — h 
program. The diagrams at each order are then evaluated using FORM [15]. The integrals 
are performed using a recursive method. The results are then combined to provide the 
needed Lagrangians at the next order. A more detailed discussion of the method and the 
underlying principles can be found in [6, 11]. 

Flavour sums in the earlier work were rather trivial to perform. Here, keeping track 
of the different terms is somewhat more tricky but all the flavour sums can be performed 
using the methods of [14]. The underlying idea is to use 

(T a AT a B) = (A) (B) - 1 (AB) , (T a A) (T a B) = (AB) - 1 (A) (B) , (6) 

for the sums over the generators T a . 

In the following, we will present our results for the coefficients of the leading logarithm 
contribution to several physical quantities. In all cases, we have the choice of expressing 
these in terms of lowest order or physical parameters, which can have quite a substantial 
effect on the convergence of the series. Following the definitions in [11] we expand a given 
observable O p h ys as 

O phys = O (1 + a x L + a 2 L 2 + ■••), (7) 
Ophys = O (l + ciL phys + c 2 Lj hys H ) , (8) 

where the chiral logarithms are defined either from the lowest-order parameters M and F 
as 

or from the physical mass M n and decay constant F n as 

M 2 ii 2 

L * E W l0g F' (10) 

7T 7T 

These are relevant for the static quantities where the mass is the only dimensionful pa- 
rameter. In general the argument of the logarithm is not uniquely determined at the level 
of leading logarithms. For the cases with more dimensionful quantities we usually use the 
more general 

r - M 2 U 2 



->M 



^St^ (ii) 



16tt 2 F2 & M 2 

where Ai is some combination of the relevant dimensionful quantities 
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i 


cti for N = 2 


ai for N = 3 


a>i for general N 


1 


-1/2 


-1/3 




2 


17/8 


27/8 


9/2N- 2 - 1/2 + 3/8 N 2 


3 


-103/24 


-3799/648 


-89/3 N- 3 + 19/3 iV" 1 - 37/24 TV - 1/12 iV 3 


4 


24367/1152 


146657/2592 


2015/8 iV~ 4 - 773/12 N' 2 + 193/18 + 121/288 N 2 
+41/72 N 4 


5 


-8821/144 


-27470059/186624 


-38684/15 N- 5 + 6633/10 N~ 3 - 59303/1080 iV" 1 
-5077/1440 N - 11327/4320 N 3 - 8743/34560 N 5 


6* 


1922964667 
6220800 


12902773163 
9331200 


7329919/240 N~ 6 - 1652293/240 iV" 4 
-4910303/15552 N~ 2 + 205365409/972000 
-69368761/7776000 N 2 + 14222209/2592000 N 4 
+3778133/3110400 N 6 



Table 1: The coefficients Oj of the leading logarithm L 1 up to % = 6 for the physical meson 
mass for the physical cases N = 2 and iV = 3 as well as for general N. 



i 


Cj for N = 2 


d for N = 3 


Cj for general TV 


1 


-1/2 


-1/3 


-N- 1 


2 


7/8 


163/72 


7/2N- 2 -3/2 + 3/8iV 2 


3 


211/48 


9329/648 


-109/6 iV~ 3 + 71/6 N- 1 - 55/24 iV + 2/3 N 3 


4 


21547/1152 


638861/7776 


3185/24 N~ 4 - 1067/12 N~ 2 + 1795/72 
-677/288 N 2 + 77/72 N 4 


5 


179341/2304 


86965829 
186624 


-150877/120 N~ 5 + 47473/60 iV~ 3 - 23407/120 N' 1 
+41713/1440 N + 4891/8640 N 3 + 57557/34560 TV 5 


6* 


2086024177 
6220800 


38355806767 
13996800 


229179/16 N~ 6 - 392659/48 N~ 4 + 106873049/77760 iV~ 2 
+2800699/60750 - 277103161/7776000 N 2 
+68361593/5184000 A^ 4 + 8001833/3110400 TV 6 



Table 2: The coefficients q of the leading logarithm L l phys up to i = 6 for the physical 
meson mass for the physical cases N = 2 and iV = 3 as well as for general N. 
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4 Mass 



The mass is known fully to one [16] and two loops [14]. We have calculated the leading 
logarithms to six-loop order here. For this we have computed the generic two-point function 
of the fields in all four parametrizations and extracted the mass as well as the wave 
function renormalization, which will be needed later. The result is expressed in the form 
of (7) and (8) with O p h ys = M% and O = M 2 . The first six coefficients a, and q of the 
expansions of the physical mass are listed in Tables 1 and 2. The coefficients for N = 2 
agree with the results from [6, 11, 12] and with the one- and two-loop results from [14, 16]. 

Since the calculation is very time consuming, the sixth order has only been checked 
with two of the four parametrizations in (4). Throughout the paper, results with this 
limitation are marked by an asterisk next to the number that labels the order. 

It is rather clear from the expressions that there is a pattern in the powers of N that 
appear. They always jump by powers of 2. Similar steps can be seen in all results quoted 
in this paper. This is due to the SU(N) group structure of all flavour traces that need to 
be evaluated as is proven in general in Appendix A. 

We can use our results to check the convergence of the two expansions. In Fig. 1 for 
N = 2 and Fig. 2 for N = 3, the input values chosen are F = 0.090 GeV for the expansion 
in terms of L and F w = 0.0922 GeV for the expansion in L p hy S as well as \i = 0.77 GeV. 
The convergence is somewhat worse for N = 3 than for N = 2. 

5 Decay constant 

The decay constant F n is defined by 

<0|jJU* 8 (p)> = iV2F^S ab (12) 

for a meson corresponding to the quark flavour combination iq r y^T a q and the axial current 
g7 At 7 5 T b g. Note that F n is equal for all mesons since we are in the equal mass limit. The 
decay constant is known fully to one [16] and two loops [14]. Here, we evaluate the leading 
logarithms to six loops. 

We need to evaluate a matrix-element with one external axial current and one incoming 
meson. The diagrams required for the wave function renormalization were already done 
in the mass calculation in the previous section. We thus need to evaluate all relevant 
one-particle-irreducible (1PI) diagrams with an external a a . Up to the more complicated 
group theory the calculation is the same as in our earlier work. 

We give the first six coefficients for both leading logarithm series with O p h ys = F n and 
Oo = F in Tables 3 and 4. Note that once F v is known as a function of F, we can express 
all observables as a function of the physical M 2 and F^. We already used this to get the 
coefficients q in Tables 2 and 4 from the corresponding aj. 

We have plotted in Figs. 3 and 4 the expansion in terms of the unrenormalized quantities 
and in terms of the physical quantities for N = 2 and N = 3 respectively. In both cases 
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Figure 1: The contribution of the leading logarithms to M 2 /M 2 order by order for F = 
0.090 GeV, F n = 0.0922 GeV, /i = 0.77 GeV and N = 2. The left panel shows the 
expansion in L keeping F fixed, the right panel the expansion in L p h ys keeping F n fixed. 
Plots similar to Fig. 1 in [12]. 
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Figure 2: The contribution of the leading logarithms to M 2 /M 2 order by order for F = 
0.090 GeV, F w = 0.0922 GeV, fi = 0.77 GeV and N = 3. The left panel shows the 
expansion in L keeping F fixed, the right panel the expansion in L p h ys keeping F n fixed. 



i 


di for N = 2 


dj for N = 3 


cii for general N 


1 


1 


3/2 


1/2N 


2 


-5/4 


-35/16 


-1/2 - 3/16 iV 2 


3 


83/24 


293/36 


23/127V- 1 + l/4iV + l/4iV 3 


4 


-3013/288 


413359 
13824 


-139/12iV- 2 + 7/54 - 523/576 N 2 - 3511/13824iV 4 


5 


2060147 
51840 


96197471 
622080 


22357/240 A^ 3 - 5063/648 iV" 1 + 16157/5184 N 
+26237/17280 N 3 + 5885/13824 iV 5 


6* 


69228787 
466560 


932532830269 
1343692800 


-41296/45 N~ A + 5690093/58320 N~ 2 
+14622631/9331200 - 1944182341/186624000 N 2 
-945730747/373248000 iV 4 - 81119291/149299200 iV 6 



Table 3: The coefficients of the leading logarithm L l up to i = 6 for the decay constant 
F n for the physical cases N = 2 and iV = 3 as well as for general N. 



i 


d for N = 2 


d for N = 3 


Cj for general N 


1 
2 
3 
4 
5 

6* 


1 

5/4 
13/12 
-577/288 
-14137/810 

37737751 
466560 


3/2 
45/16 
131/36 

113471 
13824 
65712649 
622080 

889506647989 
1343692800 


1/2N 
5/16N 2 

-1/3 N' 1 + 1/8N + 1/8 N 3 

9/4 N~ 2 + 209/576 TV 2 - 229/108 - 1639/13824 iV 4 

-1097/60 iV" 3 + 11095/648 N' 1 - 40225/5184 N 
+2137/2880 N 3 - 679/1536A^ 5 

6745/36 iV" 4 - 9274909/58320 iV~ 2 +611736991/9331200 
-3858946741/ 1866240007V 2 + 440983853/373248000 TV 4 
-127342211/149299200 N 6 



Table 4: The coefficients q of the leading logarithm L l phys up to i = 6 for the decay constant 
F n for the physical cases N = 2 and iV = 3 as well as for general N. 

we get convergence but it is better for the expansion in physical quantities. It is also much 
better for N = 2 than for N = 3. 



6 Vacuum expectation value 

The expression for the leading logarithms of the vacuum expectation value (VEV) follows 
from the definition 

Vy S = (0|-Jo|0), (13) 

where j$ is the QCD current associated with the scalar external source s introduced in (3) 
with the singlet generator normalized to 1. In terms of quarks the definition is 

{o\m\o) = ( 14 ) 
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Figure 3: The contribution of the leading logarithms to F n /F order by order for F = 
0.090 GeV, F n = 0.0922 GeV, // = 0.77 GeV and N = 2. The left panel shows the 
expansion in L keeping F fixed, the right panel the expansion in L phys keeping F n fixed. 
Plots similar to Fig. 2 in [12]. 
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Figure 4: The contribution of the leading logarithms to F n /F order by order for F = 
0.090 GeV, F n = 0.0922 GeV, \i = 0.77 GeV and N = 3. The left panel shows the 
expansion in L keeping F fixed, the right panel the expansion in L p \ iys keeping F n fixed. 
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At lowest order, Vph ys = Vq — —2BF 2 . The VEV is known fully to one [16] and two loops 
[14]. Here we evaluate the leading logarithms to seven loops. 

The first seven coefficients of the expansions in (7) and (8) for O p h ys = ^phys and Oq = Vo 
are given in Tables 5 and 6, respectively. 

We have plotted in Figures 5 and 6 the expansion in terms of the unrenormalized 
quantities and in terms of the physical quantities for N = 2 and N = 3, respectively. 
In both cases we get a good convergence but it is excellent for the expansion in physical 
quantities. 



7 Vector form factor 

We turn now to the vector form factor which is defined by 

(r(Pf)\j c v,M(p>)) = (T c {T b T a -T a T b ))(p f+Pt ),F v [(p f - Pt ) 2 } 



(15) 



The vector current is q / j fJi T c q. It is known fully in two- and three-flavour ChPT to one 
[2, 17] and two loops [18, 19]. Here we calculate the leading logarithms in the equal mass 
case to five loops. 

The procedure to find the leading logarithms is entirely the same as in the earlier work 
[11, 12] with the modifications needed for the more complicated flavour structure. We 
express the result in terms of t = t/M 2 and the logarithm (11) with a scale Ai 2 that is 
some combination of t and M 2 To fifth order we find 



F v (t) = l + L M N/12 1 



+ L 



M 



i (1/12 + N' l /\Q) +P (A^/288) 



T 3 



t (-1/12 iV- 1 + 29/72 N + 1 79/2592 N' 6 ) +P (-137/648 N - 5/2592 AT 



+ t 6 {N/32 + 5/5184 N a 

+ L 4 M t (17/18 N~ 2 - 617/576 + 1111/768 N 2 + 12011/124416 iV 4 ) 

+ t 2 (445/7776 - 79861/103680 TV 2 - 5519/414720 N 4 ) 

+ P (79459/777600 N 2 + 407/172800 N 4 ) + t 4 (-7V 2 /345600 + 7V 4 /518400) 



+ L 5 



M 



t (-391/45iV~ 3 + 134797/12960 N^ 1 - 5429753/933120 N + 5232541/1244160 N 3 

+ 14129/93312 N 5 ) + P (-181271/233280 N' 1 + 194389/405000 N 

- 262760809/116640000 TV 3 - 61724321/1866240000 N 5 ) + P (35/384iV" 1 

+ 1261489/23328000 TV + 35449669/90720000 iV 3 + 129189077/26127360000 iV 5 ) 



+ T (-226531/4320000 N - 58095211/1632960000 Ar - 935713/4354560000 
+ P (871/115200 AT + 545009/181440000 A^ 3 + 126059/6531840000 A^ 5 ) 

Note that Fy(Q) = 1 as it should be. 



(16) 
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i 


Oj for TV = 2 


di for iV = 3 


a,i for general N 


1 


3/2 


8/3 


-N~ l + N 


2 


-9/8 


-4/3 


3/2 N- 2 - 3/2 


3 


9/2 


988/81 


-20/3 N~ 3 + 22/3 A^ 1 - 7/6 iV + 1/2 iV 3 


4 


-1285/128 


-5660/243 


1025/24iV~ 4 - 205/4iV- 2 + 175/16 - 55/24iV 2 - 5/48 iV 4 


5 


46 


399563 
1944 


-350iV- 5 + 2188/5 N~ 3 - 12539/120 iV" 1 + 1321/80 N 






-373/960 N 3 + 737/960 N 5 


6* 


1305605 
9216 


242777185 
419904 


2490019/720 iV" 6 - 3137701/720 N~ 4 








+12971623/12960 N~ 2 - 1295581/12960 

+154399/51840 N 2 - 277697/69120 N 4 - 68761/207360 N 6 








7* 


153149887 
226800 


137725650367 
27556200 


-12489752/315 N~ 7 + 15424312/315 N~ 5 
-79037542/8505iV- 3 - 129606331/567000 N' 1 
+66338023/324000 N - 19464419/3402000 iV 3 
+67022189/13608000 iV 5 + 4453133/2721600 N 7 



Table 5: The coefficients of the leading logarithm L % up to i = 7 for the vacuum 
expectation value Vph ys for the physical cases N = 2 and N = 3 as well as for general N. 



i 


d for N = 2 


d for N = 3 


Cj for general N 


1 


3/2 


8/3 


—N" 1 + N 


2 


21/8 


68/9 


1/2 iV~ 2 - 3/2 + N 2 


3 


75/16 


1720/81 


-7/6 N- 3 + 7/3 N- 1 - 13/6 N + N 3 


4 


1023/128 


26881/486 


109/24 N~ 4 - 103/12 Af" 2 + 277/48 
-127/48 iV 2 + 11/12 N 4 


5 


2669/256 


82861/729 


-637/24 N- 5 + 5587/120 N~ 3 - 57887/2160 iV -1 
+9241/1080 N - 5263/2160iV 3 + 179/270 N 5 


6* 


48029 
138240 


67564919 
1049760 


150877/720 A^ 6 - 49505/144 N~ 4 + 46879/288 N~ 2 
-378373/12960 + 9427/10368 N 2 - 2741/6912 N 4 
+14701/103680 N 6 


7* 


6924628769 
87091200 


966193799261 
881798400 


-229179/112 N~ 7 + 440981/140 N~ b 
-967456169/816480AT- 3 - 986277601/163296000 iV" 1 
+3916037663/31104000 N - 56820907057/1306368000 N 3 
+7411227769/1306368000 N 5 - 19637251/26127360 N 7 



Table 6: The coefficients q of the leading logarithm L l phys up to i = 7 for the vacuum 
expectation value Vph ys for the physical cases N = 2 and N = 3 as well as for general N. 
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Figure 5: The contribution of the leading logarithms to Vph ys /Vo order by order for F = 
0.090 GeV, F n = 0.0922 GeV, fi = 0.77 GeV and N — 2. The left panel shows the 
expansion in L keeping F fixed, the right panel the expansion in L phys keeping F n fixed. 
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Figure 6: The contribution of the leading logarithms to Vph ys /Vo order by order for F = 
0.090 GeV, F n = 0.0922 GeV, fi = 0.77 GeV and N = 3. The left panel shows the 
expansion in L keeping F fixed, the right panel the expansion in L p h ys keeping F n fixed. 
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i 


d for N = 2 


d for iV = 3 


Q for general N 


1 


1 


3/2 


N/2 


2 


2 


31/8 


l/2 + 3/8iV 2 


3 


853/108 


877/48 


-1/2 N- 1 + 29/12 N + 179/432 iV 3 


4 


50513/1296 


824171/6912 


17/3 iV~ 2 - 617/96 + 1111/128 iV 2 + 12011/20736 iV 4 


5 


120401/648 


33850135/41472 


-782/15 N~ 3 + 134797/2160 iV" 1 - 5429753/155520 N 
+5232541/207360 N 3 + 14129/15552 N 5 



Table 7: The coefficients Cj of the leading logarithm £p hys up to i = 5 for the vector radius 
(r 2 )y for the physical cases N = 2 and iV = 3 as well as for general N in units of M 2 . 



i 


q for N = 2 


q for N = 3 


Cj for general iV 


1 











2 


1/72 


1/32 


iV 2 /288 


3 


-71/162 


-593/864 


-137/648 N - 5/2592 N 3 


4 


-25169/7776 


-1978981/248832 


445/7776 - 79861/103680 iV 2 - 5519/414720 iV 4 


5 


-1349303/72900 


75784701173 
1119744000 


-181271/233280 N' 1 + 194389/405000 N 
-262760809/116640000 N 3 
—61724321/ 1866240000 N 5 



Table 8: The coefficients q of the leading logarithm L l phys up to z = 5 for the curvature cy 
for the physical cases N = 2 and N = 3 as well as for general N in units of M 2 . 

The formula in the massless case is much simpler. The logarithm is now a bit more 
unique. We replace M 2 by — t and define 

Taking the limit M 2 — > (which implies F n F), we get from (16): 

F°(t) = 1 + K t (N/12) + K 2 (A^ 2 /288) + K*(N/32 + 5/5184 iV 3 ) 
+ fQ 4 (-iV7345600 + iV 4 /518400) 

+ ^(871/115200 N + 545009/181440000 N 3 + 126059/6531840000 N 5 ) . (18) 

We close this section with giving the expansion for the radius and curvature of the 
vector form factor defined by 

F v {t) = l + hr 2 ) v t + c v t 2 + -- - . (19) 
6 

The coefficients q for the expansion in physical quantities are given in Tables 7 and 8 in 
units of M 2 . The result up to two- loop order agrees with the LL extracted from the full 
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i 


d for N = 2 


d for N = 3 


Q for general iV 


1 


-1 


-2/3 


-27V- 1 


2 


31/8 


569/72 


23/2 N~ 2 - 7/2 + 9/8 N 2 


3 


65/6 


3205/81 


-260/3 iV~ 3 + 133/3 N' 1 - 95/12 iV + 23/12 N 3 


4 


76307/1152 


2330311/7776 


19801/24iV- 4 - 5731/12iV- 2 + 1033/9 
-3331/288 N 2 + 295/72 N 4 


5 


375263/1152 


65426359/31104 


-189077/20 N~ 5 + 82642/15 N~ 3 - 1352663/1080 iV" 1 
+63967/360 N - 35647/4320 iV 3 + 48487/5760 N 5 



Table 9: The coefficients q of the leading logarithm L l phys up to i = 5 for the scalar form 
factor at zero momentum transfer Fs(0) for the physical cases N = 2 and N = 3 as well 
as for general N. 

two-loop calculation [18]. We do not present numerical results for the vector form factor 
since these are dominated by large higher-order contributions, see, e.g., [2, 18]. 

All the results presented in this section agree for iV = 2 up to fifth order with the 
findings of [12]. 



8 Scalar form factor 

The (singlet) scalar form factor is defined by 

(PM\-3o\F(pi)) = FsVPf-Pi)*l 



(20) 



where again we have normalized the scalar current generator to one. It is known fully in 
two- and three-flavour ChPT to one [2, 17] and two loops [18, 20]. Here we calculate the 
leading logarithms in the equal mass case to five loops. 

As opposed to the vector case, the scalar form factor is not normalized to one, such that 
we also need to specify the LL expansion of -Fs(O). The coefficients q for the expansion 
in terms of physical logarithms for Oq = 2B and O p h ys = -^5(0) are given in Table 9. 
The momentum dependent part, Fs{t) = Fs(t)/Fs(0), can again be expressed in terms of 
t = t/M 2 and the logarithm defined in (11). To fifth order we find 



F s (t) = l + L M N/2i 



+ L 



M 



i (-1/2 - 5/18 A^) +P 43/144 iV : 



+ L 



M 



+ L 4 



M 



t{5/6N~ l -3/8N - 473/2592 N ) + r (-91/216 N - 227/1296 iV 3 ) + t 6 143/864 N 
t (-13/6 iV" 2 - 1751/1944 - 10529/5184iV 2 - 2117/3456 iV 4 ) 



r (4645/3888 + 245537/388800 N 2 + 27103/259200 N 



+ t 



J3 



-196121/388800 N 2 — 57061/388800 iV ) + t (1129/57600^ + 580837/6220800 iV 
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i 


d for N = 2 


d for iV = 3 


Cj for general N 


1 


6 


9 


3iV 


2 


-29/3 


-18 


-3-5/3iV 2 


3 


-581/54 


-1663/48 


57V- 1 - 9/4 N - 473/432 N 3 


4 


-75301/648 


-2147363/5184 


-13N~ 2 - 1751/324 - 10529/864iV 2 - 2117/576 iV 4 


5 


5482247 
9720 


535098163 
186624 


417/5 iV" 3 +327877/4860 N' 1 - 1550429/38880 N 
-28557851/777600 iV 3 - 3800759/518400 iV 5 



Table 10: The coefficients q of the leading logarithm L l phjs up to i = 5 for the the scalar 
radius (r 2 )s for the physical cases N = 2 and iV = 3 as well as for general N. 



i 


d for N = 2 


d for N = 3 


Cj for general N 


1 











2 


43/36 


43/16 


43/144 iV 2 


3 


-727/324 


-863/144 


-91/216 N - 227/1296 iV 3 


4 


4369/810 


2386939/155520 


4645/3888 + 245537/388800 N 2 + 27103/259200 iV 4 


5 


16871641 
2916000 


1130937893 
55987200 


-222149/58320 N' 1 - 858337/466560 N 
+410235883/233280000 iV 3 
-39345049/466560000 N 5 



Table 11: The coefficients q of the leading logarithm L l phys up to i = 5 for the curvature 
cs for the physical cases N = 2 and iV = 3 as well as for general N. 



+ L 5 M t (139/10 A^" 3 + 327877/29160 N' 1 - 1550429/233280 N - 28557851/4665600 N 3 

- 3800759/3110400 iV 5 ) + t 2 (-222149/58320 N' 1 - 858337/466560 N 
+ 410235883/233280000 TV 3 - 39345049/466560000 N 5 ) 
+ t 3 (324253/233280 N - 7699463/38880000 N 3 + 26029871/311040000 N 5 ) 
+ t 4 (-1129/144000 JV - 357457/1296000 iV 3 - 18692191/186624000 iV 5 ) 

+ ? (315439/25920000 iV 3 + 48727189/933120000 iV 5 ) 



(21) 



As for the vector form factor, we can also give our result for the radius and the curvature, 
which are defined as 

F s (t) = l + ];(r 2 ) s t + cst 2 + --- . (22) 
o 

The coefficients q for the expansion in physical quantities are given in Tables 10 and 11. 

All the results presented in this section agree for = 2 up to fourth order with the 
findings of [11]. 
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9 Meson- meson scattering 

The amplitude for general meson-meson scattering is defined from 

(0 c (P3)0 d (p4)out|0 a ( Pl )0 b (p 2 )in) = i(27r) 4 S\p 3 +p A - Pl -p 2 )M(s,t,u) , (23) 

with the Mandelstam variables 

s = {pi+P2) 2 , t = ( Pl -p 3 ) 2 , u ={p l -p i f. (24) 

It has been calculated at one- loop order in [21]. The two- loop calculation has been per- 
formed together with two other symmetry breaking patterns in [22]. 

The structure of the SU(N) meson-meson scattering amplitude has been derived in full 
generality in [21, 22]. It can be expressed in terms of two invariant amplitudes B(s,t,u) 
and C(s, t, u) as 

M{s,t,u) = [(T a T b T c T d ) + (T a T d T c T b )] B{s,t,u) 

+ [(T a T c T d T b ) + (T a T b T d T c )] B(t,u,s) 
+ ^T a T d T b T c ) + (T a T c T b T d )] B(u,s,t) 

+ 5 ab 5 cd C{s, t, u) + 5 ac 5 bd C{t, u, s) + 5 ad 5 bc C{u, s, t) . (25) 

The T a are the generators of SU(N) normalized as (T a T b ) = 5 ab . Crossing symmetry 
implies 

B(s,t,u) = B(u,t,s) , C{s,t,u) = C{s,u,t) . (26) 

In the case of iV = 2 the traces over four generators evaluate to products of Kronecker 
deltas such that the structure of the amplitude is reduced to the well-known expression 

M(s, t, u) = 5 ab 5 cd A{s, t, u) + 5 ac 5 bd A{t, u, s) + 5 ad 5 bc A{u, s, t) , (27) 

with 

A(s, t, u) = C(s, t, u) + B(s, t, u) + B{t, u, s) - B(u, s, t) . (28) 

This is of course the structure of the 7T7t scattering amplitude in two-flavour ChPT. 

We have calculated the LL contribution to the two invariant amplitudes to fifth order. 
In order to make the symmetries (26) explicit, we have expressed B(s,t,u) in terms of 

t = t/Ml and A su = (s - u)/M* , (29) 

and C(s,t,u) in terms of 

s = s/M% and K ut = {u - t) /Ml . (30) 
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We also express it in terms of the more general logarithm (11). Our result for general N 
reads 

F 2 

-^B(s,t,u) = l-t'/2 + L A 4-2/3(3iV- 1 - N) - 5/12 M + 1/16 NP + 1/48 NA 2 SU ] 

7T 

+ L 2 M [1/36(684N~ 2 + 12 + 29iV 2 ) - l/36(54iV~ 2 + 156 + 17N 2 )t 

+ 1/288(486 + 29iV 2 )P - 1/288(42 - lljV 2 )A 2 u - 5/1152(54 + N 2 )? 

- 5/384(6 + N 2 )tA 2 u ] 

+ L 3 M [- 1/6480(1 1880007V" 3 - 224160iV" 1 - 54388iV - 7461iV 3 ) 
+ l/12960(90720iV" 3 + 40200iV _1 - 192278iV - 7391iV 3 )t 

- l/51840(33840iV" 1 - 3179527V - 7127iV 3 )P 
+ 1/51840(2040(W _1 - 178887V + 3219iV 3 )A 2 u 

- l/51840(56769iV + 5027V 3 )? 3 - l/51840(9403iV + 1646iV 3 )iA 2 u 
+ 17/207360(711iV + 8N 3 )t A + l/207360(4122iV + 653iV 3 )?A 2 u 

- l/69120(43iV - 29iV 3 )A 4 J 

+ L 4 v( [-l/777600(-1710720000A^ 4 + 272984400iV~ 2 + 7733020 - 24562906iV 2 

- 1393823iV 4 ) + 1/3110400(-2083968007V" 4 - 216817200iV~ 2 + 32764628 

- 120616232iV 2 - 2240271iV 4 )t + l/777600(10376100iV- 2 + 11270867 

+ 12959626iV 2 + 141465iV 4 )t 2 + l/3110400(-9954000iV~ 2 + 4785940 - 1937608iV 2 
+ 310421iV 4 )A 2 u _ l/12441600(14175000iV- 2 + 90102306 + 45489100iV 2 
+ 198731iV 4 )i 3 - l/12441600(4725000iV" 2 + 4860662 + 4068564iV 2 
+ 725363iV 4 )?A 2 u + 1/24883200(33470280 + 9600665iV 2 + 520847V 4 )? 4 

- 1/2764800(281520 - 173014iV 2 - 27275iV 4 )t 2 A 2 u - 1/24883200(528120 

+ 125223iV 2 - 33181iV 4 )A 4 u - 1/24883200(2694708 + 395685iV 2 + 1270iV 4 )F 

- 1/12441600(418230 + 63783iV 2 + 7363iV 4 )?A 2 u 

- 1/8294400(135360 + 25767iV 2 + 2642iV 4 )?A 4 J 

+ L^[-l/233280000(7018928640000iV~ 5 - 894502656000iV~ 3 - 283189944960iV~ 1 

+ 72507663308iV - 22041000184iV 3 - 690252879 iV 5 ) + l/3265920000(2508879744000iV~ 5 
+ 2387314944000iV" 3 - 1788779140400iV~ 1 + 55061393444iV - 302932752254iV 3 

- 3125842333A^ 5 )t + l/1866240000(54336096000iV~ 3 - 23570392560A^~ 1 + 4485329438AT 

- 1556519119iV 3 + 304085556iV 5 )A 2 u - l/13063680000(1174116384000iV~ 3 

- 743306538800 N' 1 - 1129493851946iV - 538364942807iV 3 - 3202580448^)? 
+ 1/26127360000(40960080000A^ -3 + 14810807200iV- 1 - 46526628758iV 

- 14697445505iV 3 - 25442917527V 5 )tA 2 su + l/26127360000(122880240000iV~ 3 

- 16992114720iV' 1 - 1115282964722iV - 281556881291iV 3 - 615480176JV 5 )? 



17 



+ 1/14929920000(5298746407V~ 1 - 971377232iV - 249290704iV 3 + 44246931iV 5 )A 4 u 

+ 1/52254720000(14905099440A^ _1 - 17891241888iV + 7135232732iV 3 

+ 1114519563iV 5 )?A 2 u - l/104509440000(17859033360iV" 1 - 9646200689 12N 

- 176725206936iV 3 - 477061005iV 5 )t 4 - l/52254720000(3722683128iV + 675880271iV 3 

+ 64699403iV 5 )tA 4 u - l/26127360000(4775058204iV + 535939563iV 3 + 59878229iV 5 )t 3 A 2 u 

- l/52254720000(53926623792iV + 7406736715iV 3 + 10009775 N 5 )t 5 

- l/5971968000(1608228iV + 172954iV 3 - 42207iV 5 )A^ + l/41803776000(35608896iV 

+ 27012324iV 3 + 47416857V 5 )PA 4 U + 1/ 13934592000(86939676 N + 292918307V 3 

+ 1761661/V 5 )t 4 A 2 u + 1/209018880000(7967262600/V + 1179147888/V 3 + 20295437V 5 )? 5 ] 

(31) 



and 



1V± -K 

+ L M [2N- 2 + 3/8s 2 + l/8A 2 t ] 

+ L 2 M [-2(7N- 3 - N- 1 ) - 1/36(108./V _1 - 37iV)5 - 13/36 N§ 2 
+ 5/24iVA 2 t + 55/192 TVs" 3 + 5/192 TVs A^] 

+ L 3 V! [l/810(112860iV^ 4 + 5580iV~ 2 - 1173 + 1022iV 2 ) 

+ 1/3240(66240A^" 2 - 35174 + 18417V 2 )s" + l/2592(4104iV~ 2 + 9256 + 1661iV 2 )s" 2 
+ l/12960(6840iV" 2 + 3936 + 3811iV 2 )A 2 t - 1/12960(19044 + 3929iV 2 )s 3 
+ 1/6480(1998 + 353iV 2 )sA 2 t + 1/23040(3762 + 4177iV 2 )s 4 
+ 1/11520(1254 + 49iV 2 )s 2 A 2 t + 1/23040(418 + 13iV 2 )A 4 i ] 

+ L 4 v( [-l/25920(41990400iV~ 5 + 6609600iV~ 3 - 4545600iV _1 + 365524iV - 112385iV 3 ) 

- l/1555200(342921600iV" 3 - 85806200iV~' 1 + 26129734iV - 1244537iV 3 )s 

- l/3110400(21772800iV~ 3 - 31914800AT 1 - 112473847V - 6159957V 3 )s 2 

- l/622080(1451520iV~ 3 + 90000iV _1 - 1035376iV - 245073iV 3 )A 2 t 
+ l/6220800(12042900iV~ 1 + 1797031iV + 2858510iV 3 )s 3 

- l/2073600(864300iV~ 1 - 2339087iV - 167184iV 3 )sA 2 t 

- l/12441600(10710972iV + 2610871iV 3 )s 4 + l/103680(46775iV + 1366iV 3 )s 2 A 2 t 
+ l/2488320(197684iV +4431iV 3 )A 4 t + 1/33177600(38386447V + 3531083iV 3 )s 5 
+ l/19906560(275724iV + 4553iV 3 )sA 4 t + l/9953280(285444iV + 5093iV 3 )s 3 A 2 t ] 

+ Z/tj[l/81648000(1754648179200iV- 6 + 506658499200iV- 4 - 280670106528iV" 2 

+ 54081789480 - 2534550216iV 2 + 999897691iV 4 ) + l/1632960000(4317509952000iV" 4 

- 1246702968000iV~ 2 + 33240243520 - 42063222570iV 2 + 131586137V 4 )S 

+ l/6531840000(147311136000iV" 4 + 93672190080iV~ 2 - 67579015880 + 41982479569iV 2 
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+ 34358375297V 4 )A 2 4 + l/6531840000(441933408000iV~ 4 - 5249942342407V -2 

+ 414501298952 + 85767984257iV 2 + 6718221609iV 4 )s 2 + l/4354560000(21832020000iV~ 2 

- 10349834500 + 12661599885iV 2 + 429983342iV 4 )sA 2 t 

- l/13063680000(238082755680iV- 2 + 224361657012 + 458283032997V 2 

+ 836998070iV 4 )s 3 + l/26127360000(2833329240iV" 2 + 1543460180 + 59757444987V 2 
+ 97700803iV 4 )A 4 t + 1/26127360000(16999975440A^" 2 + 11020206360 + 32781452788iV 2 
+ 783206303iV 4 )s 2 A 2 4 + l/13063680000(12749981580iV- 2 + 117034928874 
+ 16745354989iV 2 + 4551556991iV 4 )s 4 + 1/10450944000(1480020120 + 1322536120iV 2 
+ 14756647iV 4 )s 3 A 2 t + 1/20901888000(1604443896 + 1525971760iV 2 + 18544639iV 4 )sA 4 t 

- 1/34836480000(37975069080 + 25358850464iV 2 + 4923041771iV 4 )s 5 

+ 1/580608000(20634624 + 4848765iV 2 + 66668iV 4 )s 4 A 2 t + 1/2322432000(77354172 

+ 15213300iV 2 + 129899A^ 4 )s 2 A 4 < + 1/52254720000(118623420 + 27848862iV 2 

+ 216847iV 4 )A^ + 1/34836480000(2583793620 + 2952657684iV 2 + 2122427329iV 4 )s 6 ] . 



The five-loop contribution is of the same calculational complexity as the mass to the sixth 
order and as the latter, has only been checked in two of the four parametrizations in (4). 
Up to fourth order and for N = 2, the amplitude agrees with A(s,t,u) from [11]. 

It is well known that for 577(2) tttt scattering, the amplitude can be decomposed into 
three amplitudes corresponding to intermediate states of fixed isospin 0, 1, or 2. This 
decomposition can be generalized to arbitrary values of N, where one finds seven different 
intermediate states. The corresponding amplitudes are obtained from the above invariant 
amplitudes as [22] 



+ C{t,u,s) + C(u,s,t) , 
T A = N[-B(s, t, u) + B{t, u, s)} + C{t, u, s) - C{u, s, t) , 
T S a = C(t, u, s) - C(u, s, t) , 
T A s = C(t, it, s) - C(u, s, t) , 
T ss = 2B(u, s, t) + C(t, u, s) + C{u, s, t) , 

T AA = - 2B(u, s, t) + C(t, u, s) + C{u, s, t) . (33) 



(32) 




+ (iV 2 - l)C(s, t, u) + C(t, u, s) + C(u, s, t) , 
[N - 4 J [B(s, t, u) + B(t, u, s)] - ^-B(u, s, t) 




The full scattering amplitude is then built up from these as 



M(s,t,u) = '%2Tj{s,t,u)Pj, 



(34) 
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i 


a for N = 2 


a for iV = 3 


aft for N = 3 


aft* for N = 2 


aft ti for iV = 3 


1 


9/2 


358/51 


59/21 


-3/2 


-14/9 


2 


857/42 


28487/612 


3505/252 


-31/6 


-955/108 


3 


153211/1512 


7143269/22032 


751735/9072 


-7103/216 


-255265/3888 


4 


41581/84 


98674513/44064 


26921179/51840 


-7802/45 


-6097649/12960 


5* 


139816697 
56700 


165016031929 
10575360 


228658804229 
65318400 


3326573 
3375 


83190853 
23040 



Table 12: The coefficients Cj of the leading logarithm L l phys up to i = 5 for the s-wave 
scattering lengths for the physical cases N = 2 and iV = 3. 



i 


af for N = 2 


af for N = 3 


af* = af a for iV = 3 


1 


2 


13/6 


2/3 


2 


791/36 


941/36 


50/9 


3 


8528/81 


665171/3888 


25481/972 


4 


2291903/3888 


678064381/559872 


51822143/279936 


5* 


894986647/291600 


234732737339/27993600 


3480221279/2799360 



Table 13: The coefficients Cj of the leading logarithm £p hys up to i = 5 for the p-wave 
scattering lengths for the physical cases N = 2 and iV = 3. 

where Pj are the respective projection operators. Since their explicit form is rather lengthy, 
we do not reproduce it here and refer the interested reader to [22]. For N = 2, the channels 
with J = S, SA, AS, AA do not exist and for N = 3 the channel with J = AA is not present. 
Each channel can be projected on partial waves by 

i r 1 

T e J (s) = -^J ^d{cos6)P e {cos6)Tj{s,t,u). (35) 

Expanding these around threshold in powers of q 2 = s/4 — M 2 leads to the definition of 
the threshold parameters: 

Re r/(a) = q 2 \aj + q 2 bj + g% J + ■••), (36) 

where aj are the scattering lengths and bj the slope parameters. We have calculated all the 
s- and p-wave scattering lengths. Note that for each channel, only one of the two partial 
waves is non-zero. For N = 2, only three channels contribute and the scattering length are 
more commonly denoted by 

<2q = Oq , a\ = af , a l = a o S ■ (37) 

For N = 3, T AA is the only channel that vanishes, such that there are six scattering lengths. 
The tree-level expressions for general N as well as for the physical cases N = 2 and N = 3 
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Figure 7: The contribution of the leading logarithms to a^/a 1 ^^ order by order for = 
0.0922 GeV, /x = 0.77 GeV and N = 3. The left panel shows the expansion in L p h ys keeping 
Fn fixed. The right panel shows the first two orders together with with the equivalent full 
results of [22], labeled NLO and NNLO. 



are given by [21, 22]. 



I, tree 
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1 
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Ml 
16ttF2 ' 



(3* 



The scattering lengths from channels that do not contribute for N = 2 or N = 3 have 
been omitted. af A,tree and af 4 ' 1 ^ 66 vanish at tree level, but the higher-order contributions 
are non-zero. 



As usual, the LL expansion can be written in the form of (8), where now 0$ 



J,tree 



There is, however, an exception: since the tree-level contribution to af A and a AS vanishes, 
the series is written in these cases as 



SA 



AS 



M, 



lQnFl 



^(Cl-Lp hyS + C 2 Lp hys H ) . 



(39) 
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The corresponding coefficients for the two physical cases are listed in Table 12 for the 
s-wave and in Table 13 for the p-wave scattering lengths. 

Our results for the scattering lengths agree up to two- loops with [22]. Furthermore, for 
N = 2 the LL contributions to and a^ s have been calculated up to fourth order in [11] 
and are in agreement with the present results. 

We show the convergence of the leading logarithm part of the scattering length for 
N = 3 in the left panel of Fig. 7. On the right side, we show the one- and two-loop LL 
compared with the full calculation of [22]. One can see that the leading logarithms are 
about half of the full correction. 

10 77 — > tttt and pion polarizabilities 

The process 77 — > inr has been calculated in ChPT to one loop in [23, 24]. Already before 
the p % Lagrangian was explicitly known, 77 — > tt°tt° has been worked out to two loops in [25] 
and 77 — > 7r + 7r~ in [26, 27]. Both calculations were redone and the 0(p 6 ) counter terms 
added explicitly in [28, 29]. While the leading contribution to the charged process comes 
from tree-level diagrams, the neutral process only starts at the one-loop level. In that case, 
knowing at least the leading logarithms at higher orders is a particularly welcome check 
for the convergence of the chiral expansion. Starting at the four-loop order, there is also 
a doubly anomalous contribution to the amplitude. However, it only affects sub-leading 
logarithms and is therefore not relevant in the present context. 

We have calculated the leading logarithms for the amplitude with two different vectors 
to two different mesons in general but the expressions are extremely lengthy. We thus 
restrict ourselves to the simpler case where both vectors are coupling to the current q r ) lx T°q 
and denote them as r y c . For the numerical results we treat only the case where the vectors 
correspond to photons and are on-shell. 

The 7 C 7 C — y 7i a 7i b scattering amplitude is defined from the matrix element 

(7r a (p 1 )7r b (p 2 )out|7 c (fc 1 )7 c (fc 2 )in) = z(2tt) + p 2 - h - k 2 )T abc (s, t, u) , (40) 

with 

s = {pi+p 2 ) 2 , t={p 1 -k 1 ) 2 , u =( Pl -k 2 ) 2 , (41) 

and 

T abc (s,t,u) = e 2 ete^(s,t,u) . (42) 

The polarization vectors for the external vector are ei,e 2 and we have added an overall 
coupling constant e to the vectors. We will consider the process for both vectors off-shell 
such that the amplitude also depends on k\ and k\. Since both vectors carry the same 
flavour index, the amplitude must satisfy 

KV, v = k»V^ = Q. (43) 
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This follows from the Ward identity. If the vectors had different flavour there could have 
been an equal time part in the Ward identity but it vanishes here. As a result the amplitude 
can be be decomposed into gauge invariant quantities as 



V» v {s, t, u) = A(s, t, u)Tjf" + B(s, t, u)T£ u + C{s, t, u)T£ v + D(s, t, u)T£ v + E(s, t, u)T^ 



pv 
5 ' 

(44) 



where 



Tipv — fa ■ k 2 g liu — k 2[1 fa v , 

T 2fiU = fa ■ k 2 A fl A u + ^(t- u)k 2fi A l/ - ^(t - u)A^fa v - ^-(t - u) 2 g^ u , 
1)ifiv — fa ■ k 2 fa^k 2u k-yk 2jl k 2v k 2 fa^,fa v + k^k 2 g^ v , 
T 4 ^u = fa ■ k 2 k llA A u - fa ■ k 2 A^k 2v - k\k 2pb A v + fa^A^fay 

+ \( t ~ + k l)d^ ~\^~ u ) k ^ k ^ ~ 2^ ~~ u )fa»fav > 
T 5flu = k\klA^A u - i(t - iijklA^v + ^(t - u)klk lfM A v - -(t - uffa^fay , (45) 

with A = pi — p 2 . For k\ = k\ = 0, T 5fMU is equivalent to T 3fMU , such that only four 
independent quantities remain, which are identical to those given in [25] up to normalization 
factors. For on-shell photons, one has in addition t\ ■ fa = e 2 ■ k 2 = 0, which only leaves 
Tx^v and T 2flv . The T iiliV satisfy the identities 

kiTi^y = k 2 Tip, v = , (46) 

which can be readily checked using 

-A- fa = A- k 2 = \{t-u) . (47) 

The flavour structure of the amplitude consists of all possible traces of T a ,T b and T c . 
But the total amplitude is symmetric under (pi, a) 4-> (p2,b) and under (fa, ei) -H- (k 2 ,e 2 ). 
Using charge conjugation 1 one can prove that the amplitude must be separately invariant 
under (a <->■ b). 

In contrast to all the earlier quantities discussed in this paper, we must here consider 
several classes of diagrams in order to calculate the physical 77 — > tttt amplitude. There are 
one-particle-reducible contributions here other than those taken care of by wave function 
renormalization. The three needed types of diagrams are depicted in Figure 8. On the one 
hand, there is the direct contribution involving the 777^ leading-order vertex, on the other 
hand, there are two types of diagrams involving the 77T7T vertex twice. Loop contributions 
to the latter vertex have been calculated already for the vector form factor and the results 



1 This is valid at least for the real part of the amplitude. 
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Figure 8: The three types of diagrams that contribute to 77 — > nn. The gray circles stand 
for the sum over all one-particle-irreducible diagrams. Note that the intermediate meson 
propagator in the two rightmost diagrams is off-shell. 



from there can be reused. The corrections to the two-point function are incorporated by 
using the physical propagator for off-shell momenta, which is given by 



1 



l + £(p 2 )p 2 -M 2 



with S(p 2 ) 



s( P 2 ) - e(m 2 ; 

p 2 _ M 2 



(46 



where £(p 2 ) denotes the self-energy of the meson. The contributions from the diagrams 
involving a propagator can be written in a compact form because these terms depend 
explicitly on the vector form factor. In particular, because of the LSZ theorem, the residue 
of the amplitude when the intermediate propagator is on-shell must contain the vector 
form factor twice. 

We introduce a notation for the flavour traces that is a little shorter to write: 



^ {rparpbrpcrpc\ _j_ {rpbrparpcrpc\ 

t 3 =(T a T b )(T c T c ) , 



t 2 = (T*T C T b T C ) ^ 

U =(T a T c )(T b T c ) . 



(49) 



That the traces in t\ only appear together is due to the symmetry under a <-> b. 

A large part of the amplitude is contained in the generalized Born amplitude because 
of the argument given above, with 



e^V^ n =(t 1 -2t 2 )F v (kl)F v (kl) 



2e 1 ■ e 2 



-A ■ e 1 A • e 2 + A ■ ei&4 ■ e 2 - k 2 ■ e x A ■ e 2 + k 2 ■ 6\ k x ■ e 2 

t ~ M 2 _ 
-A • ei A ■ e 2 - A • eifci ■ e 2 + k 2 ■ e x A • e 2 + k 2 ■ e x k x ■ e 2 

u-Ml 



(50) 
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The remainder of the amplitude now has no poles in the t or u channel. The pole amplitude 
also has the commutator structure expected from tree-level couplings to external vectors, 
visible in T abc . The expression for Fy(t) in terms of the LLs can be found in (16). The 
generalized Born amplitude can be decomposed in the functions defined in (44), but is 
somewhat simpler in the form given in (50). 
We write the full amplitude now as 

V,u =V* orn + A/ Ml T 1(IV + B/M* TV + C/M* T^ v + D/M* T 4 , u + E/M* T^ v . (51) 

The factors of M 2 are introduced to make the functions dimensionless. The partial ampli- 
tudes we write as functions of 

h = kl/M*, k 2 = k 2 2 /M 2 , k 12 = h-k 2 /M 2 n , A tu =(t-u)/M 2 . (52) 

Each of the amplitudes we then write as 

A = A^ Ll, + A^L 3 M + A&L A M + ■■■. (53) 

The leading logarithms at one-loop order are already fully contained in (50). The two- loop 
leading logarithms are quite simple: 

B {2) = (1/72 N 2 + l/12)ti - 1/6 1 2 + 1/9 Nt 3 - 1/36 Nt 4 , 

A m = (k 12 -2)B^ 2 \ C® = -B<?\ D® = E® = 0. (54) 

The third-order expressions are still reasonable in full generality: 

A® = t 4 [l/216 (30kf 2 + 9k 2 k 12 - hk\ + 9kJ: l2 - 2k x k 2 - 5k 2 + 64k 12 + 54k 2 + 54&i - 104) 

+ 1/1296(30^ + Uk 2 k 12 ~H + 14jfcijfei2 + 4^2 - kf - 32k 12 + 50k 2 + 50h + 88)N 2 ] 
+ t 3 [1/432(30^2 + 9k 2 k 12 - hk\ + 9k x k 12 - 2k x k 2 - hk\ + 64fci 2 + 54fc 2 + 54&i - 104) 
+ 1/2592(30A;?2 - 3\k 2 k 12 - 31k\ - 2,lk x k l2 - 50hk 2 - 37k\ + 548fci 2 - 1072)iV 2 ] 
+ t 2 [1/1296(210^2 + 21^2^12 - 299A| + 21fcifcia + 58%^ - 299A; 2 + 608fci 2 
+ 2042^2 + 2042fci + 1832)iV - 1/1296(2^12 + 8fcf + 2k x k l2 - k x k 2 + 8kj - 26k l2 
- 52k 2 - 52fci + 12)iV 3 ] 

+ ti [-1/1296(60^ - 3*2*12 - 142fc| - 3hk 12 + 32kik 2 - U2kj + 208&i 2 + 940fc 2 
+ 940&i + 1072)JV - 1/5184(30^ 2 + 15* 2 *i2 - 13*1 + 15*i*i 2 + Uk^ - 13*? - 40*i 2 
+ 160^ + 160fci + 184)iV 3 ] , 
S (3) = t 4 [l/216(30jfei2 + 29fc 2 + 29A;i - 20) + l/648(15*i 2 + 13* 2 + 13Jfei - 58)iV 2 ] 

+ t 3 [l/432(30*i 2 + 29fc 2 + 29*i - 20) + l/2592(30fei 2 + 53* 2 + 53&i + 464)iV 2 ] 

+ h [1/I296(210fci2 + 293fc 2 + 293&1 - 1996)iV + l/432(* 2 + *i - W)N 3 ] 

+ ti [-1/I296(60fci2 + 103* 2 + 103*i - 968) N - l/5184(30*i 2 + 29* 2 + 29*i - 236)iV 3 ] 
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C (3) = t 4 [-l/216(34jfei2 + 17fc 2 + 17fci - 36) - l/324(3fci 2 + fe 2 + ^1 - 27)iV 2 ] 

+ t 3 [- 1/432(34^12 + 17fc 2 + 17fci - 44) - l/2592(210&i 2 + 121fc 2 + 121fci + 452)iV 2 ] 

+ t 2 [- 1/1296(106^12 + 77k 2 + 77fci - 1484)iV - l/1296(fc 2 + h - 18)iV 3 ] 

+ *i [1/I296(2fci2 + 13^2 + 13fei - 676)iV - l/5184(10fci 2 + 7fc 2 + 7fci + 200)iV 3 ] 

L>( 3 ) = t 4 A ta (5/216 + 5/1296 iV 2 ) + t 3 A te (5/432 + 5/2592 iV 2 ) + i 2 A ta (53/1296 N + 1/2592 iV 3 ) 
+ tiAt u (-19/1296iV - 1/864 N 3 ) 
= t 4 (7/54 + 13/648 N 2 ) + t 3 (7/108 + 11/648 iV 2 ) + i 2 (23/162iV) 

+ h(-25/648N - 1/216 iV 3 ) . (55) 

The fourth-order expression is very long. We therefore only quote the on-shell case with 
k\ = _k% — €i • ki — e 2 • k 2 — 0, where the amplitude is reduced to the contributions from A 
and B: 

= t 4 [-85/324(fci 2 - 2)N~ 1 + l/777600(129fci 2 A 2 tu - 40908fc? 2 + 12292A t 2 u - 573518fc 2 2 
+ 3289918fei2 - 2678500)iV - l/1555200(174£;i 2 A 2 U + 14452fc? 2 + 1027A 2 U 

- 79808fc 2 2 - 67142^12 - 75300)iV 3 ] 

+ t 3 [43/324(^12 - 2)N~ 1 + l/1555200(5772fci 2 A 2 U + 346056fc? 2 + 14281A 2 M - 187274£; 2 2 
+ 1943624fei2 - 2697000)iV + l/3110400(453A;i2 A 2 U + 37544fc? 2 - 806A 2 n + 6624fc 2 2 
+ 1042176A:i2 - 1985600)iV 3 ] 

+ * 2 [-19/54(jfei2 - 2)N~ 2 - 1/388800(1881^12 A 2 tu + 128988fc? 2 + 663A 2 M + 128748fc 2 2 

- 444098^12 - 862700) - l/1555200(4824fei 2 A 2 tu + 326352fc? 2 + 50627A 2 U - 1077658£: 2 2 

- 5276492^12 - 3115200)iV 2 - l/62208(39A 2 n + 280fc 2 2 - 5810&12 + 6228)iV 4 ] 

+ *i [19/108(^2 - 2)iV" 2 + 1/777600(1881 A;i 2 A 2 n + 128988fc? 2 + 663A 2 U + 128748fe 2 2 

- 460898^2 - 829100) + l/518400(1139fci 2 A 2 tu + 69072fc 3 2 + 10597A 2 M - 253738fc 2 2 

- 478512fci 2 - 925400)iV 2 + 1/6220800(219^2 A 2 tu + 12512fc? 2 + 5012A 2 M - 174348fc 2 2 
+ 48248^12 - 20800)iV 4 ] 

fi (4) = U [~ 85/324 N- 1 + 1/777600(129A 2 U - 40908fe 2 2 + 564866^i 2 - 362950)iV 

- 1/777600(87A 2 U + 7226fe 2 2 - 83402fei 2 + 170475)iV 3 ] 

+ t 3 [43/324 N- 1 + 1/777600(2886A 2 M + 173028fe 2 2 - 323231&12 + 216900)iV 

+ 1/3110400(453A 2 U + 37544fc 2 2 - 65888fci 2 + 812200)iV 3 ] 

+ t 2 [- 19/54 N~ 2 - 1/388800(1881 A 2 tu + 128988fe 2 2 - 403776^2 + 68650) 

- l/777600(2412A t 2 u + 163176fe 2 2 - 1391027fei 2 + 4833750)iV 2 + 29/10368(2fei 2 - 29)iV 4 ] 
+ h [19/108 N~ 2 + 1/777600(1881 A 2 tu + 128988fc? 2 - 403776fci 2 + 51850) 

+ 1/1555200(3417A 2 M + 207216fc 2 2 - 1310482fci 2 + 4768800)iV 2 

+ 1/6220800(219A 2 U + 12512fe 2 2 - 162724fci 2 + 691400)iV 4 ] (56) 
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The two- loop leading logarithms 2 for N = 2 agree with those of [28, 29] for A and B. 
Note that the in the formulas given there also contain a logarithm denoted by £. 

We can now use these results to find the polarizabilities. These are defined from the 
helicity amplitudes, 

H++ ~Ml + ~MT B ^ H+ ~ Mis B ' (5?) 



at fixed t = M\ 



2. 



jfH +T (s, t = M 2 n ) = ( ai ± ft) + -(a 2 ± fo) + 0(s 2 ) . (5* 



The leading terms are the dipole, the next-to-leading terms the quadrupole polarizabilities. 
They can then be expanded as 

a < ±A =16^( C - + S + °< M ->)' (M) 
At one-loop order the leading logarithms vanish. This was shown for iV = 2, 3 in the earlier 
works [23, 24, 30]. The two- loop LLs for N = 2 can be extracted most easily from the 
expressions in [28, 29]. The conclusion is that the only terms containing LLs are the d\ + 
and from [28, 29] we get 

di + ^°) = I logV /Ml) , d 1+ (n+) = i logV/M*) • (60) 
Alternatively we can write the expression for the polarizabilities using our notation as 

7T 

t 

«2±/3 2 =^^c^ hys . (61) 

7r i 

The resulting coefficients are given in Table 14. Note that they should not be confused 
with the coefficients q± in (59). We only quote the results for N = 2. The general iV 
results depend on how one extends the charge matrix to more flavours. 

In Table 15 we have given the numerical LL contributions together with the full two- loop 
results in units of 10 -4 fm 3 for the dipole and 10 -4 fm 3 for the quadrupole polarizabilities. 
The results are somewhat mixed. For the phenomenologically most relevant case, — Pi, 
the LL are always small and well below the uncertainty quoted in [28, 29]. Also a\ + Pi for 
the neutral pion is not much affected but the other + /3, obtain significant corrections 
to their actual value. The estimate in [28] of the p 8 contributions from omega exchange to 
«2 + P2 neutral is —0.25 in the same units, again much larger than the LL. The charged 
«j + Pi obtain larger relative corrections since the two-loop estimate of [29] is very small. 
Note that the chiral logarithm contributions in [28, 29] also include nonleading logarithms 
which is why those numbers differ from the ones in Table 15. 



2 Our calculation is with a charge matrix with vanishing trace. However, the singlet part does not appear 
in the lowest-order Lagrangian, hence we get the correct result for N = 2 using Q = diag(l/2, — 1/2) rather 
than Q = diag(2/3, -1/3). 
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-8/9 
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-352/27 
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4631 
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1045 


114829 
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18 
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81 


972 


54 


2025 
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Table 14: The coefficients q of the leading logarithms for the pion polarizabilities as defined 
in (61) for the physical case N = 2. 
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0.11 











0.23 











3 


0.011 


-0.022 


0.11 


0.26 


0.10 


0.039 


-0.056 


-0.65 


4 


0.001 


-0.011 


0.017 


0.14 


0.037 


0.024 


-0.14 


-0.42 


[28, 29] 


1.1 


-1.9 


0.04 


37.6 


0.16 


5.7 


-0.001 


16.2 



Table 15: The numerical contribution of the leading logarithms at each order to the pion 
polarizabilities with the physical charged pion mass M n = 139.57 MeV, F n = 92.2 MeV 
and \x = 0.77 GeV. The units are 10 -4 fm 3 for the dipole and 10 -4 fm 5 for the quadrupole 
polarizabilities. 

11 Conclusions 

In this work we extended the earlier work on leading logarithms in effective theories and 
especially massive nonlinear sigma models to the case of SU(N) x SU(N)/SU(N) with 
equal meson masses. We presented results for the leading logarithms for up to 7 loops for 
the mass, decay constant, vacuum expectation value, vector form factor, and scalar form 
factor, as well as for a number of quantities connected with meson-meson scattering and 
77 — > tttt. When applicable we have provided results for both physically relevant cases 
N = 2 and N = 3 and for general N. In all cases we find reasonable convergence for 
N = 3, while it is even better for N = 2. 

All the results presented here have been checked through the use of several parame- 
terizations. Furthermore, we have compared to existing one- and two-loop results. The 
higher orders could be checked for iV = 2 by comparison with our previous work on the 
0(3) model. 

We have not done a general study of how well the large iV limit works. But looking at 
the coefficients in the various tables, one notes that the coefficients of the subleading terms 
are often larger than the leading coefficients even though usually not by much. Some cases 
with substantial corrections to the large iV limit can however be found. In conclusion, the 
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large iV limit is not typically a good approximation to the full result but is significantly 
better than was found for the 0{n) case, both due to the size of the coefficients and the 
fact that the suppression is now in powers of iV 2 . 

One of the motivations behind this work was the hope that knowing many of the 
leading coefficients in N would allow for an educated guess at the all-order series. We 
did not succeed in this. Our work can serve as a starting point for future studies in this 
direction. 
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A Powers of N in the results 

The formulas and tables in the main text show a clear step of 2 in the powers of the number 
of flavours N that show up. In this appendix we show that this must be the case to all 
orders. The same is true for the colour factors in QCD with iV c -colours due to gluonic 
contributions but we are not aware of a simple published proof even though it might be 
well known in the QCD perturbation theory community. 

If we had a theory with the symmetry breaking pattern U(N) x U(N)/U(N) we would 
have N 2 Goldstone bosons and the proof could be done using 't Hooft's double line notation 
with the only difference that the lines indicate flavour of the quarks rather than colour. 
The argument in [31] in the purely gluonic theory with surfaces with handles goes through 
in the same way and one concludes that only positive powers of N should appear and that 
adding a handle changes the power by two. 

Since in our case we do not have N 2 particles but rather N 2 — 1, we cannot simply take 
over the arguments from [31]. Those still determine the highest power of N at each loop 
order, though. 

In the main text we determined the leading logarithms from one-loop diagrams only. 
But the leading logarithms are in principle also determined from the diagrams with largest 
number of loops at any order. We use the fact that in these diagrams, all vertices come from 
the single trace lowest-order Lagrangian. The proof immediately generalizes if external 
fields are included. We give the proof for one-particle-irreducible diagrams only but the 
generalization should be fairly clear. 
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For a given one-particle-irreducible diagram, the number of loops Nl, vertices Ny, and 
propagators N P satisfy the relation 

Nl = Np — Ny + 1 . (62) 

For the parametrization U\ in (4), exactly one flavour trace appears for each vertex, such 
that the total number of traces in the diagram is Ny. For each of the Np propagators 
there is a sum over a flavour index, which can be evaluated using the relations 

(T a A) (T a B) = (AB) - ±{A){B), (T a AT a B) = (A)(B) - i(AB) . (63) 

We denote the total number of traces in a term by N tr and the power of N by N^. The sum 
of the two is abbreviated by N tr N = N tr + N^. Each time the first one of the relations (63) 
is used, N tT N is reduced by one. The first term of the second relation, on the other hand, 
adds one to N tr N, while the second term again subtracts one. Positive powers of iV can be 
generated by (1) = N. 

After all propagator flavour traces have been removed, a diagram with Ny vertices can 
contain terms with N tr N = Ny — Np, Ny — Np + 2, . . .. The minimal value is achieved if 
N^n has been decreased by one Np times. If the first term of the second relation in (63) has 
entered exactly once, one gets instead N tr N = Ny — (Np — 1) + 1 and so on. From (62) then 
follows that N trN is odd (even) for even (odd) Nl- The tree- level one-particle- irreducible 
diagrams clearly satisfy this since they contain no powers of A" and one flavour trace. 

For a given number of loops, the lowest number that can occur is N tr N = 1 — Nl 
using (62). The highest requires a little more work because not all contractions can increase 
NtrN- In a one-particle-irreducible loop diagram with Ny vertices, we need to use the first 
relation in (63), which only lowers N tr N, at least Ny — 1 times since at least that many 
operations have the T a in different traces. So the maximum is 

NrN = Ny- (Ny - 1) + (N P — (Ny — 1)) = Np — Ny + 2=1 + Nl . (64) 

This coincides with the maximum power derived with the double line method of [31]. 
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